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Abstract 

We show theoretically that thermomechanical effects in dye-doped nematic liquid crystals when 
illuminated by laser beams, can become important and lead to molecular reorientation at intensities 
substantially lower than that needed for optical Freedericksz transition. We propose a ID model 
that assumes homogenous intensity distribution in the plane of the layer and is capable to describe 
such a thermally induced threshold lowering. We consider a particular geometry, with a linearly 
polarized light incident perpendicularly on a layer of homeotropically aligned dye-doped nematics. 

PACS numbers: 05.45.-a, 42.70.Df, 42.65. Sf 
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I. INTRODUCTION 



Optically induced temperature changes in liquid crystals are at the origin of interesting 
nonlinear behaviors [1]. Indeed, due to the light absorption, thermal effects can change the 
physical properties of the liquid crystal, which in turn affects the light propagation in the 
medium. Among the examples of thermal effects, we can distinguish between a direct change 
of the refractive indices, often referred to as thermal indexing, and the variation with the 
temperature of other physical parameters, such as the elastic constants, which may influence 
the light propagation as well. In the latter case, thermal effects might be responsible for 
director reorientation, and are thus referred to as thermomechanical effects. 

Previously, light-induced thermomechanical effects have been largely investigated in 
cholesteric liquid crystals, where these effects were originally related with the absence of 
the right-left symmetry [2j. Then it became clear that such effects might also exist in sys- 
tems which possesses this symmetry, such as the nematic liquid crystals (NLC). As was 



shown in 
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4j, thermomechanical effects do give a contribution to the director, heat and 
Navier-Stokes equations. These thermomechanical terms are given by nonlinear combina- 
tions with respect to temperature, director and velocity gradients and represent nonlinear 
cross couplings between them. The first question which arose is to suggest an experiment 
for measuring the magnitude of thermomechanical coefficients. Indeed, it is hardly possible 
to find a setup which leads to the contribution of one of the coefficients only. Since for 
uniformly aligned nematics no thermomechanical effects are expected, the hybrid-oriented 
nematics were used in the experiments a a. There, the hydrodynamic flow appeared as 
a result of the applied temperature gradient, which allowed to find the magnitude for the 
thermomechanical coefficients. Laser induced thermomechanical effects in dye-doped ne- 
matics have been envisaged in preliminary experiments [7|. However, it was not possible 
to derive a definitive conclusion due to the main difficulty of separating thermomechanical 
contributions from the light-induced molecular torque. Another problem when performing 
these experiments is that of avoiding to approach the nematic-isotropic transition, where 
an enhancement of the nonlinear optical response of dye-doped nematics could take place 
because of other effects, such as the weakening of the anchoring jjj. 

Here, we study theoretically the thermomechanical effects which occur in uniformly 
aligned dye-doped nematic liquid crystals. This happens when the sample is illuminated 
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by a laser beam with a wavelength in the absorption band of the dye, which causes a sig- 
nificant heating of the liquid crystal layer. Note that the absorbtion is negligible for pure 
nematics, so that thermomechanical effects are significant only for dye-doped nematics. In- 
deed, in the presence of dye-doping thermal heating leads to an additional torque which acts 
onto the director together with the light-induced torque. This additional torque will help 
either to destabilize or stabilize the initial orientation which is determined by the ratio of the 
thermomechanical coefficients. We show that the additional torque might be strong enough 
to induce the molecular reorientation well before the onset of the light-induced molecular 
reorientation, the so called optical Freedericksz transition (OFT) for pure liquid crystals 



9[] and Janossy effect for dye-doped liquid crystals [10|, llll]. We account for the possible 
decrease of the reorientation threshold by a ID model that assumes homogenous intensity 
distribution in the plane of the layer (the plane wave approximation) and includes the light 
absorption into the hydrodynamic equations for the nematics. 

II. THEORETICAL MODEL 

In dye-doped nematics, contrary to the pure ones, the OFT lowering might happen not 



only because of the well known Janossy effect [10|, lllf but also because thermomechanical 
effects may add a significant contribution. Indeed, as a consequence of light absorption, 
light propagation in dye-doped nematics causes significant heating of the LC. When the 
intensity of a beam is sufficiently large and temperature gradient becomes nonzero, an 
additional force acting on the director and additional terms in the stress tensor in Navier- 
Stokes equation for the velocity appear due to thermomechanical effect {3, 4\. The equation 
for the velocity is coupled with the director equation, so any dynamical process that leads 
to director reorientation will also induce flow even in the absence of pressure gradients. It 
should be noted that the terms which describe the thermomechanical effects appear not only 
in the director but also in Navier-Stokes equations. Thus, neglecting the velocity equations 
and considering only the director equation may lead to misleading results, even though this 
simplification has been done for many studies in the context of light induced instabilities. 

In our case, we considered the full problem with the velocity equations coupled to the 
director reorientation equation. However, in order to calculate the shift of the OFT thresh- 
old due to the thermomechanical effect we employed another frequently used simplification, 
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FIG. 1: Geometry of the setup: a linearly polarized light along the x-direction is at a normal 
incidence on a nematic LC layer; no || z is the director of the initial unperturbed alignment 
(homeotropic state). 

namely that all variables depend only on one coordinate, which is transversal to the plane 
of the nematic layer (ID assumption). We thus considered a linearly-polarized plane wave 
incident perpendicularly on a layer of a dye-doped nematic of thickness L that is sand- 
wiched between two substrates of thickness d (see Fig. [1]). The cell has initially homeotropic 
alignment (with strong homeotropic anchoring at the boundaries) and is placed in a ther- 
mostage with a temperature To on both sides. The light is polarized along the x-direction 
and propagates along the positive z-axis. 

A. Heat equation 

For simplicity we also assumed that i) the attenuation of the light inside the nematic 
is small [i.e. I(z) ~ I Q , where Jo is the incident intensity]; ii) the transversal heat flow 
occurring in the (x, y) plane is neglected, i.e. the temperature profile depends only on z, 
T = T(z). Since we deal with the plane- wave approximation the following ID steady-state 
heat conductivity equations in the nematic and substrates can be written (note that an 
absorption takes place only inside the nematic): 

K s d 2 z T = 0, -d < z < or L < z < L + d (1) 
K\\d?T = —a±Io, < z < L , 

where k\\ is the parallel component of the heat conductivity tensor of the nematic, k s is the 
heat conductivity of the substrates and a± is the absorption coefficient for the ordinary light. 
We then write the boundary conditions given by continuity conditions of the temperature 
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and the heat flow at the substrate- nematic interfaces (z = and z = L): 

T s i\ z =o = Tjy\ z=0 , K s d z T s i\ z= o = K\\d z T]y\ z=0 (2) 
Tn\ z =l = T s2 \ z= l, K\\d z T N \ z=L = K s d z T s2 \ z= L- (3) 

In this simplest model the temperature profile is linear inside the substrates and has a 
parabolic form inside the nematic which is symmetric with respect to the center of the 
layer z = L/2. The maximal temperature is at the center of the layer and the maximal 
temperature difference inside the nematic is AT max = a±IoL 2 /(8k\\). The temperature 
gradient inside the nematic can be written as 

d z T N = (3p{L- 2z), with (3 = (4) 
where p = Iq/If is the incident intensity normalized to the threshold intensity of the OFT 
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7T 2 c(e± + e a )K 3 e a + C 

If = -pi p= , V = • (5) 

L 2 e a Jelr} e a 



for a dye- doped nematic 



Here e a = ey — e±_ is the dielectric anisotropy and e±_ (en) is the dielectric permittivity (at 
optical frequency) perpendicular (parallel) to n, ( phenomenologically describes the effect 
of certain dye dopants (( = in a pure LC), K3 is the bend elastic constant of the nematic 
and c is the velocity of light in the vacuum. 

It should be noted that the obtained solution (0J for the temperature profile inside the 
nematics is much simpler than in reality. The more complicated solution for Gaussian 
incident beams has been derived in [l3| from 2D heat equation which includes transversal 
dependence. It turned out that for large Gaussian beams i.e. when the spot size w is much 
larger than the thickness of the layer L, w ^> L, the maximum temperature rise, AT max , 
becomes proportional to the spot size. Such a behavior is not predicted by the ID model and 
is owing to the transversal heat flow occurring in the plane of the layer. Nevertheless, the 
realization of the plane wave approximation in the experiment is difficult but not impossible 
task. One of the way to proceed is to enlarge the Gaussian beam to a size which is much 
larger than the medium working area (which is much stronger condition then w ^> L), so 
that the intensity can be considered uniform in the central part. Another possibility is that 
of using the so-called flat-top beams [14]. Both techniques are quite easily accessible and 
allow to attain the threshold for OFT. Preliminary experiments are running at present in 
our laboratory and will be reported elsewhere. 



B. Linearized Navier-Stokes equation 

The Navier-Stokes equation for the velocity v can be written as [11 

p m (d t + v VH = -V, (p 5 i:i + TTij + T™ sc + T™) , (6) 
where p m and p are the density and the pressure of the LC, respectively, tt^ is the Ericksen 



stress tensor 



151 ] . The viscous stress tensor T^ sc in Eq. is written in terms of the 
six Leslie coefficients a, [jjj] and the thermomechanical tensor Tf- M is introduced in ji], 4|. 
Then, the incompressibility condition (p m is constant) V ■ v = and the no-slip boundary 
conditions v| 2= o,l = immediately ensure that the z component of the velocity vanishes 
v = {v x {z, t), v y (z, t), 0), so v is parallel to the plane of the layer. Moreover, all convective 
derivatives v ■ V vanish. The Navier-Stokes equation has been simplified then due to the 
fact that the director relaxation time 



differs by many orders of magnitude with the momentum diffusion time r visc = p m L 2 /'~f 1 , 
where 71 = «3 — a 2 is the rotational viscosity. Typically T V i SC ~ 10~ 6 s and r ~ Is, so the 
slow variable of the system is the evolution of the director which enslaves the flow motion 
and, thus, the inertial terms in the Eq. ([6]) can be neglected. Taking into account that the 
light is polarized in the x-direction, we need an equation for the x component only. The 
linearization of this equation around the homeotropic state (n x — 0, v x — 0) yields 

v 

{a 5 - a 2 + a 4 )^^ + a 2 d t n x - (8) 
a 10 d z n x d z T N = C(t), 

where C(t) is a function that does not depend on z and will be fixed by the boundary 
conditions. Here the first two terms on the left-hand side of Eq. ([8]) come from the viscous 
stress tensor whereas the last one is the contribution from the thermomechanical stress 
tensor with the thermomechanical coefficient aio- (Note that aio in [J] is related to the & 
from 3] as a w = (£s — ^4)/4-) The unknown function C(t) can be determined by integrating 
Eq. ([8]) across the layer. Finally, the velocity gradient d z v x can be expressed in terms of the 
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director and the temperature gradient [which is given by Eq. (j3J)] as: 

2 

d z v x = ■ x (9) 

«5 — «2 + a 4 

[C(t) - a 2 d t n x + a 10 /3p(L - 2z) d z n x ] 

C(t) = ^rla 2 J d t n x dz - 2a w f3p J n x dz\ . (10) 

V / 

C. Linearized director equation. Adiabatic elimination of the flow field 

The equation for the director n is 



71 (d t + v • V - ux) n = (7 2 ^n + h g 1M ) , (11) 

where 72 = 03 + a 2 - h is the molecular field obtained from the variational derivatives of 
the free energy density F, which consists of the elastic and the electrical parts [ijj]. The 
projection operator 8^ = 5ij — riiUj in Eq. (fTTj) ensures conservation of the normalization 
n 2 = 1. Here Aij is the symmetric strain-rate tensor and the vector N gives the rate of 
change of the director relative to the fluid. In Eq. ffTT]) g™ is the contribution of the 



thermomechanical effect to the force acting on the director 

We then linearized Eq. (|TT|) around the homeotropic state and obtained the following 
equation for n x 

lid t n x + a 2 d z v x = (12) 



d 2 z n x + ( 



7T\ 2 

pn x 



- y d z n x ■ d z T N . 



Note that the first term in the square brackets on the right-hand side of Eq. f|T2l) stems from 
the linearization of h whereas the last one after the linearization of g T . 

From here on we will use normalized time t — > t/r [where r is the director relaxation 
time, see Eq. (JTJ)], length z — > irz/L (the same symbols will be kept). We will also intro- 
duce dimensionless viscosity coefficients a i = ojj/71. Eliminating the velocity gradient from 
Eq. f|T2|) with the help of Eq. ([9]) and using the expression for the temperature gradient (@J, 
the following equation for n x can be derived 

d 2 z n x + p n x - (1 - b)d t n x - dip(-n - 2z)d z n x - (13) 

d t n x dz ^ / n x dz } = 

a 2 








where b, d\ and d 2 are dimensionless parameters denned as follows: 

b = — — ^ r , di = Slip, d 2 = 5 2 ip. 



(14) 



a 5 — a 2 + a 4 

Here ip depends on the absorption and 8\^ 2 are the algebraic combinations of the thermome 
chanical coefficients given by 



^ _ fir _ a±c(e± + e a ) 
71 2K||e av /eI?7 



— — + — , ^2 = ai . 
ct 9 2 



D. Linear stability analysis of the homeotropic state 

We look for solutions of Eq. Q13D of the form 

n x (z,t) = n x (z)e ct , 
where cr is the growth rate and obtain from Eqs. ( jl~3l) 

+ [p - cr(l - 6)] - dip(7r - 2z)d z n 2 
n x dz = . 



- I <t 

7T 



2d 2 p 



On 
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(16) 



(17) 



It should be noted that Eq. (lT7j) reduces to the classical linearized equation for the OFT when 
the flow and the thermomechanical effect are neglected by putting 6 = and d\ — d 2 — 0. 
Taking into account the boundary conditions n x \ z=Qt7T = 0, Eq. ffTTj) is solved by 
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where \Fi is the confluent hypergeometric function 16j and A is some constant which de- 
pends on the director itself 



A 



b ( 2d 2 p 
— Iff — 

7T \ a 2 



77,,, dz. 



(19) 



Substituting the solution for n x (|T8|) into Eq. (fT9]) the equation for the growth rate a versus 
incident intensity p can be derived (note that A will be cancelled). Substituting o = into 



8 



this equation, the following transcendental equation for the critical intensity p c has been 
obtained 



7T 
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FIG. 2: The critical intensity p c versus parameter 8i under assumptions that 82 = 8\ (solid line), 
82 = (dot-dashed line), 82 = —8\ (dashed line). p c = 1 corresponds to the critical intensity for 
the OFT (no thermomechanical effect). 

In the calculations, we took the material parameters for the nematic E7: K 3 = 
15.97 x 10~ 12 N , n e = 1.746, n = 1.522 (extraordinary and ordinary refractive indices), 
A = 514 nm, a± = 88 cm' 1 (absorption coefficient), Ku = 10~ 3 W/°C cm (the heat con- 
ductivity for the nematic). The calculations were made for a layer of 75 pm thickness. For 
these parameters and for the threshold value Ip = 32 W/cm 2 observed in our preliminary 



experiment, £ turns out to be ( ~ 52 which is a typical value for the enhancement factor [17 ] 
[see also Eq. fl5}]. The maximal temperature difference inside the layer is AT max ~ 20 K. 
In the calculations we used b = 0.8 and a 2 = —1.058. The thermomechanical coefficients 
and, hence, Si and S 2 are_unknown. Thus we took the typical order of their magnitudes as 
1(T 12 N/°C reported in 
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FIG. 3: Contour lines for the surface jO c (^ij^2)- The values of p c are shown on contour lines. 
(p c = 1 for 5\ = 62 = corresponds to the case without thermomechanical effect.) 

In Fig. [2] the results of numerical solution of Eq. fl20|) [or alternatively the eigenvalue 
problem (fl7j) ] are shown for three different cases when i) 62 — Si; ii) 62 = 0; iii) 82 = —Si. 
{82 = corresponds to the thermomechanical single-constant approximation.) One can see 
from this figure that p c might be indeed several times lower than that for the OFT. This 
fact is clearly demonstrated in Fig. [3] where the contour lines for the surface p c (^i,^2) are 
plotted. One is now forced to conclude that the thermomechanical effect leads indeed to the 
substantial change of the critical intensity. 

E. Thermomechanical effect due to temperature difference at the boundaries 

We have assumed so far that the temperature on both the bounding plates is the same. 
As a next step, we have analyzed the influence of the nonzero temperature difference AT 
maintained at the boundaries to the instability threshold. To obtain this effect in pure form 
we assumed that the thermomechanical effect is due to AT only and the absorption inside 
nematic is neglected. For this simple situation the temperature gradient inside nematic is 
constant. Following the similar procedure described in previous section when linearizing 
basic equations around the homeotropic state [see Eq. (|T7|) ]. the following ODE for n x has 
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been derived 



d 2 z n x + [p-a(l- b)] n x - d 3 ATd z n 



X 



(21) 



7T 



bo 
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n x dz = 



TT 







where 



d 3 = 8x1], r] 



(k s L + 2K\\d)7rK 3 ' 



(22) 



and the other quantities have the same meaning as before. After substituting o = in 
Eq. (|2~TT) the following simple formula for the critical intensity was found 



As is seen p c is always higher than the threshold for OFT and depends quadratically on 
AT. This effect is, however, small because the thermomechanical coefficient enters quadrat- 
ically to the expression for p c as well. (The order of magnitude of d 3 is 10~ 2 °C~ l for the 
parameters used in calculations.) 

III. CONCLUSIONS 

In conclusion, we have shown theoretically that thermomechanical effects might be at the 
origin of significant lowering of the OFT threshold expected for dye-doped nematic liquid 
crystals. To explain this, we have developed a simple model, assuming that all physical 
quantities depend only on the coordinate across the layer. We linearize both the director 
and Navier-Stokes equations around the basic state, to assess the change of the primary 
instability due to thermomechanical effects. The temperature gradient across the layer, 
which is induced by light itself due to absorption of the dye dopants, was calculated from 
the ID heat equation. Using a typical value for the thermomechanical coefficients, we have 
found that the effect of OFT's lowering might be explained by thermomechanical effects. We 
have also analyzed a situation when the thermomechanical effects are due to the temperature 
difference maintained at the boundaries. It turned out that in this case they always lead to 
an increase of the OFT threshold. 



p c = 1 + d 3 



a (AT) 2 



(23) 
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